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Energy,  exergy,  economic  and  environmental  (4E)  analysis  and  optimization  of  a  hybrid  solid  oxide  fuel 
cell  and  micro  gas  turbine  (SOFC-MGT)  system  for  use  as  combined  generation  of  heat  and  power  (CHP) 
is  investigated  in  this  paper.  The  hybrid  system  is  modeled  and  performance  related  results  are  validated 
using  available  data  in  literature.  Then  a  multi-objective  optimization  approach  based  on  genetic  algo¬ 
rithm  is  incorporated.  Eight  system  design  parameters  are  selected  for  the  optimization  procedure. 
System  exergy  efficiency  and  total  cost  rate  (including  capital  or  investment  cost,  operational  cost  and 
penalty  cost  of  environmental  emissions)  are  the  two  objectives.  The  effects  of  fuel  unit  cost,  capital 
investment  and  system  power  output  on  optimum  design  parameters  are  also  investigated.  It  is  observed 
that  the  most  sensitive  and  important  design  parameter  in  the  hybrid  system  is  fuel  cell  current  density 
which  has  a  significant  effect  on  the  balance  between  system  cost  and  efficiency.  The  selected  design 
point  from  the  Pareto  distribution  of  optimization  results  indicates  a  total  system  exergy  efficiency  of 
60.7%,  with  estimated  electrical  energy  cost  0.057  $kW  1  h  and  payback  period  of  about  6.3  years  for 
the  investment. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  world’s  increasing  energy  demand  and  the  environmental 
impacts  of  fossil  fuel  utilization  in  conventional  power  plants  have 
attracted  much  interest  to  the  development  of  clean  and  high  ef¬ 
ficiency  energy  systems.  Still  majorly  relying  on  fossil  fuels,  fuel 
cells  stand  as  the  most  promising  technology  for  their  ability  to 
avoid  combustion  emissions  due  to  the  fact  that  they  directly 
convert  the  chemical  energy  of  fuel  into  electricity.  High  efficiency, 
clean  operation  and  fuel  flexibility  of  fuel  cells  has  made  them  a 
very  attractive  choice  [  1  ]  for  energy  specialists.  The  solid  oxide  fuel 
cell  has  advantages  such  as  increasing  power  output  and  efficiency 
over  other  types  of  fuel  cells  in  power  generation  industries  when 
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they  are  used  as  a  bottoming  cycle  with  a  gas  turbine.  The  resulted 
hybrid  system  can  reach  the  first  law  efficiency  levels  more  than 
60%  and  overall  thermal  efficiency  more  than  80%  [2],  This  system  is 
applicable  for  small  scale  plants  by  coupling  with  micro  gas  tur¬ 
bines  (MGTs)  for  distributed  generation.  Modeling  and  optimizing 
the  hybrid  SOFC-micro  gas  turbine  system  is  a  subject  of  interest. 
Lubelli  and  Massardo  [3]  performed  a  comprehensive  study  on 
different  layouts  of  hybrid  SOFC— GT  systems  for  the  atmospheric 
and  pressurized  configurations  and  studied  their  performance 
under  the  effects  of  different  parameters.  Campanari  [4]  and  Cos- 
tamagna  [5]  predicted  the  performance  of  a  hybrid  MGT-SOFC 
system  in  both  full  load  and  part  load  conditions  for  mode  and 
speed  control.  Chan  et  al.  [6—10]  proposed  an  accurate  model  for 
the  electrochemical  behavior  of  the  fuel  cell  and  studied  the  effects 
of  various  parameters  on  the  simple  and  hybrid  SOFC  systems  in  full 
and  part  loads  which  provided  a  means  for  sizing  the  system. 
Massardo  and  Magistri  [11]  performed  an  exergy  analysis  on  the 
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Nomenclature 

W 

power  (kW) 

x,y,z 

chemical  reaction  rate  (kmol  s-1) 

ACPR 

air  compressor  pressure  ratio 

At 

heat  transfer  area  (m2) 

Subscripts 

C 

cost  rate  ($s-1) 

AC 

air  compressor 

C 

cost ($) 

AIP 

air  injection  pipe 

C 

specific  cost  ($lcg  1 ),  ($lcW  1  hr1) 

An 

anode 

Cp 

specific  heat  (kj  mol-1) 

AUX 

auxiliary 

E 

exergy  rate  (kW) 

act 

activation  voltage  loss 

e 

specific  exergy  (kj  kg-1) 

c 

cold 

F 

Faraday’s  Constant 

Ca 

cathode 

FCPR 

fuel  compressor  pressure  ratio 

Cap 

capital 

Gibbs  function  of  formation  (kj  kg-1) 

CC 

combustion  chamber 

H 

enthalpy  rate  (kW) 

Cone 

concentration  voltage  loss 

h 

enthalpy  (kj  kg-1) 

El 

electrochemical 

hf 

enthalpy  of  formation  (kj  kg-1) 

env 

environmental 

I 

electric  current  (A) 

ex 

exergetic 

i 

current  density,  interest  rate  (Am-2),  (%) 

/ 

fuel 

Kp 

chemical  equilibrium  constant 

FC 

fuel  compressor 

LHV 

lower  heating  value  (kj  kg-1) 

FPH 

fuel  pre-heater 

LMTD 

logarithmic  mean  temperature  difference  (K) 

h 

hot 

m 

mass  flow  rate  (kgs-1 ) 

HR 

heat  recovery 

N 

number  of  working  hours  per  year  (hours) 

O&M 

operation  and  maintenance 

n 

system  lifetime  (years) 

Ohm 

Ohmic  voltage  loss 

ft 

molar  flow  rate  (kmol  s-1) 

react 

reaction 

P 

pressure  (bar) 

REC 

recuperator 

P 

payback  period  (years) 

Ref 

reforming  reaction 

4 

heat  transfer  rate  (kW) 

Rev 

reversible 

R 

universal  gas  constant  (kj  kmol-1  K-1) 

Shift 

shifting  reaction 

SIC 

steam  to  carbon  ratio 

T 

turbine 

T 

temperature  (K) 

TPR 

turbine  pressure  ratio 

Greek  symbols 

U 

overall  heat  transfer  coefficient  (Wm-2  K) 

£ 

effectiveness 

Uf 

fuel  utilization  factor 

V 

efficiency,  voltage  loss 

V 

voltage,  volume  (Volt),  (m3) 

<p 

maintenance  factor 

hybrid  system  and  carried  out  a  thermo-economic  study  on  this 
system  by  proposing  cost  models.  Bavarsad  [12]  investigated  the 
hybrid  system  from  the  second  law  viewpoint  as  well  as  effects  of 
various  parameters  on  the  exergy  destruction  in  system  compo¬ 
nents.  Calise  et  al.  [13—16]  modeled  and  analyzed  the  hybrid  sys¬ 
tem  from  energy,  exergy  and  economic  aspects  and  performed  a 
single  level  thermo-economic  optimization  analysis.  Autissier  et  al. 
[17]  used  an  optimization  approach  to  design  the  hybrid  system  for 
two  objectives  of  efficiency  and  cost.  Duan  et  al.  [18]  performed  a 
parametric  optimization  of  the  hybrid  SOFC-MGT  system. 

However,  in  the  most  of  above  research  works  the  first  law  ef¬ 
ficiency  was  applied  as  the  criteria  for  evaluating  the  system  per¬ 
formance,  this  paper  proposes  four  E  (energy,  exergy,  economic, 
environmental)  analysis  and  multi-objective  optimization  for  esti¬ 
mating  the  values  of  interested  design  parameters  (decision  vari¬ 
ables)  of  a  hybrid  SOFC  and  micro  gas  turbine  CHP  system.  Two 
objective  functions  used  in  optimization  procedure  were  the  exergy 
efficiency  and  system  total  cost  rate.  The  environmental  penalty 
cost  caused  by  CO2  emissions  was  also  added  to  the  total  cost.  An 
optimization  procedure  based  on  genetic  algorithm  has  been 
incorporated  to  reach  the  optimum  design  parameters  with  regard 
to  a  set  of  constraints.  The  results  were  discussed  and  an  optimum 
solution  point  was  selected  from  the  Pareto  front.  The  exergy 
destruction  in  system  components  and  effects  of  optimization  on 
the  exergy  destruction  has  been  studied,  in  a  next  step  the  sensi¬ 
tivity  analysis  of  the  optimum  solutions  to  different  parameters 


were  performed.  The  followings  are  the  contribution  of  this  paper 
into  the  subject: 

•  Simultaneous  energy,  exergy,  economic  and  environmental 
(4E)  modeling  and  analysis  of  the  hybrid  SOFC-MGT  CHP 
system  was  performed  to  predict  the  system  performance,  for 
generating  combined  heat  and  power. 

•  The  optimal  values  of  system  design  parameters  were  esti¬ 
mated  using  new  objective  functions,  new  design  parameters 
and  a  list  of  constraints. 

•  Performing  sensitivity  analysis  of  the  system  to  study  the 
change  in  optimum  values  of  design  parameters  with  variation 
of  fuel  cost  as  well  as  the  capital  investment  and  system  power 
output. 


2.  Plant  description 

The  SOFC-MGT  system  layout  is  depicted  in  Fig.  1.  This  config¬ 
uration  includes  all  critical  thermal  processes  of  a  direct  coupled 
internal  reforming  SOFC-MGT  hybrid  system  in  a  simplified 
manner  and  is  close  to  the  first  real  constructed  system  of  its  kind 
[2,4]. 

The  system  consists  of  the  following  processes: 

Air  (point  1 )  and  CH4  fuel  (point  9)  at  known  inlet  conditions 
pass  through  compressors  (AC  and  FC).  These  streams  (points  2  and 
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Desulfurized 
Fuel  9 


AC:  Air  Compressor 
FC:  Fuel  Compressor 
REC:  Recuperator 
FPH:  Fuel  Pre-Heater 
AIP:  Air  injection  Pipe 
G:  Generator 
GT:  Gas  Turbine 
CC:  Combustion  chamber 
PR:  Pre-Reformer 
IRSOFC:  Internal  Reforming 
Solid  Oxide  Fuel  Cell 
HR:  Heat  Recovery 


10)  are  preheated  in  recuperator  (REC)  and  fuel  preheater  (FPH) 
respectively  using  the  recovered  energy  from  turbine  exhaust 
gases.  Since  only  hydrogen  can  participate  in  the  electrochemical 
reaction  of  the  fuel  cells,  natural  gas  should  be  reformed  to 
hydrogen  rich  products  in  presence  of  steam  at  high  temperature 
and  pressures.  By  design,  the  internal  structure  of  SOFC  allows  for  a 
large  portion  of  fuel  conversion  to  take  place  at  the  pressure  and 
temperature  of  the  fuel  cell  itself  by  means  of  internal  reforming 
capability.  The  heat  needed  to  drive  the  reforming  process  can  be 
produced  by  exothermic  electrochemical  reaction.  It  is  also  vital  in 
direct  reforming  that  some  percentage  of  inlet  methane  be  con¬ 
verted  to  hydrogen  before  entering  the  anode  compartment  in  a 
pre-reforming  reactor.  The  steam  needed  for  these  reactions  can  be 
provided  by  recirculation  of  the  anode  outlet  products  to  the 
entering  stream.  By  this  way  the  need  to  a  costly  separate  steam 
generator  would  be  eliminated.  On  the  other  side,  air  after  recu¬ 
peration  (at  point  3)  passes  through  another  pre-heater  (Air  in¬ 
jection  pipe,  AIP)  and  then  is  sent  to  the  cathode.  Then  the 
electrochemical  reaction  takes  place  on  the  three  phase  boundaries 
of  the  electrodes  where  a  large  amount  of  electrical  and  thermal 
energy  is  released.  In  this  type  of  fuel  cells,  the  air  before  reaching 
the  electrode  surface  is  heated  again  in  the  air  injection  pipes  by  the 
heat  of  combustion  reaction.  This  occurs  when  the  outlet  streams 
from  the  anode  and  excess  air  from  the  cathode  get  mixed  (Fig.  2) 
where  not  reacted  fuel  and  all  combustible  products  are  fully 
oxidized.  For  this  reason  the  air  injection  pipes  are  modeled  by  use 
of  a  virtual  heat  exchanger.  Hot  gases  after  this  point  (5)  are 
expanded  trough  a  radial  gas  turbine  where  additional  power  is 
extracted.  The  combustion  products  which  one  produced  in  the 
combustion  compartment  (CC)  and  passed  through  the  turbine  are 
mainly  sent  to  the  recuperator  (about  95%  in  stream  7)  while  the 
rest  about  5  percents  is  sent  to  fuel  pre-heater  (stream  8).  The 
exhaust  flow  (stream  14)  still  has  sufficient  thermal  energy  to  be 
recovered  for  domestic  water  heating  consumptions  in  a  heat  re¬ 
covery  exchanger  (HR).  The  desulfurized  inlet  fuel  permits  the 
system  exhaust  gases  to  be  cooled  down  to  as  low  as  50  °C  without 
the  risk  of  formation  of  sulfuric  acid  and  stack  corrosion  [2], 


3.  Thermal  modeling 

For  thermal  modeling  of  the  system  in  steady  state  condition 
the  average  values  of  the  thermodynamic  parameters  at  each 
component  were  applied.  In  this  model  the  gradients  of  tem¬ 
perature  and  pressure  in  the  fuel  cell  as  well  as  other  components 
were  neglected  and  the  equilibrium  outlet  temperature  of  fuel 
cell  stack  was  considered  as  its  working  temperature.  For  thermo¬ 
physical  properties  of  gases,  a  temperature  dependent  specific 
heat  model  based  on  empirical  polynomials  for  ideal  gas  was 
applied. 

3.1.  Compressors  and  gas  turbine 


The  outlet  temperatures  of  fuel  and  air  compressors  as  well  as 
the  gas  turbine  were  computed  from  the  following  equations  by  try 
and  error  scheme. 


Fig.  2.  Internal  reforming  tubular  SOFC  schematic. 
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For  air  and  fuel  compressors: 


/c'-‘"-^in<Ac  pr)  n) 

T, 

T,  o 

/  Cp,fudf  =  ^ln(FCPR)  (2) 

Tg 

For  turbine: 

t6 

J  =  ^rln(TPR)  (3) 

t5 

Cp,gas,avgln  (J)  =  R%ln(TPR)  (4) 

3.2.  Heat  exchangers 


used  to  determine  the  temperature  and  thermodynamic  state  of 
exiting  flows  from  the  stack. 

3.3.1.  Chemical  equilibrium  model 

It  is  assumed  that  the  following  reactions  take  place  inside  the 
stack  and  each  individual  cell.  From  which  at  the  anode  compart¬ 
ment,  the  reforming  and  shifting  reactions  were  supposed  to  be  in 
chemical  equilibrium.  In  order  to  determine  the  composition  of 
outlet  gases  from  the  cell,  the  following  reactions  were  studied: 

xCH4  +  xH2  0  <-*■  xCO  +  3xH2  Reforming  (8.1) 

yCO  +  y H2  0  <->yC 02  +  yH2  Shifting  ( 8.2 ) 

zH2  +  ZQ  02)  <->zH20  Electrochemical  (8.3) 

where:  x,  y  and  z,  denote  the  rate  of  reforming,  shifting  and  overall 
electrochemical  reactions  respectively. 

The  constants  of  equilibrium  for  each  reaction  can  be  described 
in  terms  of  partial  pressures  of  reactants  and  products  or  concen¬ 
tration  of  these  species  as  follows: 


The  design  specifications  of  recuperator  and  fuel  pre-heater 
(plate  fin  type  heat  exchanger)  were  obtained  by  using  e-NTU 
method  for  cross  flow  configuration  without  mixing  of  hot  and  cold 
fluids  [19],  For  known  design  effectiveness,  the  unknown  temper¬ 
atures  could  be  computed  from  the  following  equations  in  which 
the  specific  heat  values  were  estimated  at  the  hot  (h)  or  cold  (c)  side 
temperatures. 


Kp,R  = 


PcqP3h2 

PctuP^o' 


Kp,s  = 


PcQ2PH2 

PC0PH20 


(9) 


The  values  of  equilibrium  constants  can  also  be  obtained  from  the 
experimental  data  in  terms  of  temperature  dependent  polynomials 
as  follows,  for  which  the  constants  can  be  found  in  Ref.  [7], 


„  _  Ch(Thl  -  Th2)  CciTq  -  Tq) 
*  ■  CminCT/n  -  Tcl)  Cmin(Thl  -  Tci ) 


(5) 


The  heat  transfer  surface  area  for  each  exchanger  was  then  ob¬ 
tained  from  the  definition  of  number  of  transfer  units  (NTU)  and  its 
relation  with  effectiveness  (e).  Appropriate  equations  and  data  in 
Refs.  [20-24]  were  used  to  compute  the  overall  heat  transfer 
coefficient. 

The  air  injection  pipes  were  also  modeled  as  a  cross  flow  heat 
exchanger  and  were  assumed  to  be  an  integrated  part  of  the  stack. 
Furthermore  a  shell  and  tube  type  heat  exchanger  was  selected  for 
heat  recovery,  for  which  the  following  relations  were  used  to 
determine  the  heat  transfer  surface  area: 


At 


Q 

UMmean 


(6) 


ATm  =  0.9ATLMtd  (7) 

The  average  value  of  overall  heat  transfer  coefficient  was  consid¬ 
ered  to  be  100  Wm-2  IT1  [23], 


3.3.  SOFC  internal  reforming 

The  solid  oxide  fuel  cell  was  based  on  a  tubular,  internal 
reforming  design,  for  which  the  performance  related  data  and 
design  parameters  were  taken  from  [6,13,15].  The  fuel  cell  modeling 
was  performed  in  three  steps.  At  the  first  step,  all  internal  chemical 
reactions  and  their  rates  were  determined  by  the  chemical  equi¬ 
librium  model  by  which  the  kinetic  conditions  and  molar  flow  rates 
of  species  were  found.  The  output  cell  voltage  and  electrical  power 
was  then  computed  by  means  of  an  electrochemical  model 
considering  all  major  voltage  losses  in  the  cell  based  on  its  working 
parameters.  Finally,  a  conventional  first  law  energy  balance  was 


log(/Cp)  =  AT4  +  BT3  +CT2  +DT  +  E  (1 0) 

Due  to  the  high  rate  of  reforming  reaction  in  the  forward  di¬ 
rection  and  the  significant  value  of  its  equilibrium  constant,  this 
reaction  can  be  assumed  to  be  driven  to  completion.  This 
assumption  avoids  the  formation  of  a  highly  nonlinear  system  of 
equations  [13], 

Having  known  the  rate  values  and  establishing  the  equilibrium 
relations  in  terms  of  molar  concentrations  taken  form  kinetics  of 
each  reaction  and  by  knowing  that  the  molar  rate  of  hydrogen 
participated  at  the  electrochemical  reaction  was  controlled  by  the 
fuel  utilization  factor  (which  is  a  design  parameter),  a  system  of 
three  equations  with  three  unknowns  (x,  y,  z)  was  formed  as  the 
following: 


"01.1, in 

(11.1) 

(nH2,in  +  3x + y  -  z)  (hCo2,in  +  y) 

(%2o,in  -  X  -  y  +  z)  (hco  in  +  x  -  y) 

(11.2) 

Uf(3x+y) 

(11.3) 

By  solving  this  system  of  equations,  the  molar  composition  of 
exiting  flows  was  determined.  Moreover,  by  the  known  rate  of 
electrochemical  reaction,  we  could  estimate  the  produced  electrical 
current  of  the  stack.  Where  F  denotes  the  Faraday’s  constant. 

I  =  2  Fz  (12) 

3.3.2.  Electrochemical  model 

The  ideal  reversible  voltage  of  the  cell  can  be  computed  from  the 
Nemst  Eq.  (13)  at  its  working  temperature  and  pressure  of  re¬ 
actants.  Where  Ag°  denotes  the  change  of  Gibbs  function  for  the 
electrochemical  reaction  at  standard  thermodynamic  conditions. 


298 


S.  Sanaye,  A. . 


bi  /  Journal  of  Power  Sources  247  (2014)  294—306 


(13) 


The  actual  output  voltage  was  then  found  by  subtracting  the 
voltage  drops  due  to  ohmic  losses,  activation  over  potential  and 
concentration  loss  from  the  ideal  cell  voltage  as  follows: 

^  =  Wev  —  Vohm  ~  (VdctAn  +  Vact,Ca)  ~  Vconc  (14) 

These  losses  are  calculated  based  on  cell  working  parameters 
and  its  design  specifications  using  the  relations  and  data  in  Ref.  [6]. 


3.3.3.  The  energy  balance 

The  first  law  of  thermodynamics  and  energy  balance  for  the 
whole  stack  was  established  as  follows  for  the  enthalpy  rate  of 
input  and  output  streams,  the  electrical  power  produced  and  the 
total  heat  rate  released  due  to  internal  reactions: 


Hin  +  Hreact  =  W  +  Hout 


(15.1) 


Hreact  =:-.(zAhe|  +  yA/lshift  +  xAhRef) 


(15.2) 


Hin  =  Eh,in/  CpjdT 


system.  In  this  study,  the  whole  model  was  developed  as  a  com¬ 
puter  code  in  MATLAB  software.  Since  there  was  a  high  integrity 
and  interdependency  between  the  system  components,  direct 
calculation  of  many  parameters  was  not  possible  at  the  first  step. 
Some  parameters  such  as  inlet  temperature  of  gases  into  the  heat 
exchangers  were  not  known  at  this  stage  and  they  were  estimated 
through  a  try  and  error  procedure.  As  the  main  strategy  of  modeling 
algorithm,  a  guessed  value  was  assigned  to  the  SOFC  working 
temperature  (I4).  Based  on  this  temperature  and  other  known 
design  parameters,  as  well  as  an  estimated  value  for  pressure  los¬ 
ses,  the  constants  of  equilibrium  in  Eq.  (10)  for  reactions  in  SOFC 
modeling  were  determined.  Then  with  the  known  assumed  inlet 
flow  rates  of  air  and  fuel,  the  system  of  Eq.  (9)  was  solved  to  find  the 
rates  of  internal  reactions  (Eq.  (8)).  Based  on  these  known  rates,  the 
exiting  gases  molar  compositions  were  estimated  from  the  kinetics 
relations  (Eq.  (8)).  Then  the  combustion  relations  (Eq.  (16))  were 
applied  and  the  molar  composition  and  temperature  of  turbine 
inlet  stream  were  determined.  These  data  provided  all  system  flow 
temperature  values.  In  order  to  evaluate  the  guessed  SOFC  working 
temperature  (T4)  value,  the  energy  balance  equation  (Eq.  (15))  was 
applied.  Then  the  SOFC  outlet  temperature  (T4)  was  modified  by  the 
energy  balance  relations  until  reaching  a  convergence  value  for  the 
SOFC  temperature  (T4).  The  flowchart  in  Fig.  3  shows  the  whole 
process  of  system  thermal  modeling.  Due  to  need  for  large  number 
of  iterations  in  computation,  the  system  modeling  was  very  time 
consuming.  Therefore  the  modified  weight  coefficients  were 
applied  to  accelerate  the  program  convergence. 


Hout  =  E"i,out  /  CpjdT  (15.4) 


W  =  VI  (15.5) 

3.4.  Combustion  compartment 

It  is  assumed  that  the  following  reactions  take  place  in  cells  after 
mixing  exiting  gases  from  the  anode  and  exiting  excess  air  from  the 
cathode.  By  these  reactions,  it  was  assumed  that  all  the  combustible 
products  exiting  the  cell  anode  are  oxidized  completely.  The  molar 
composition  of  combustion  products  was  defined  easily  through 
the  reaction  kinetics.  To  find  the  gas  mixture  temperature,  an  en¬ 
ergy  balance  for  adiabatic  combustion  chamber  was  established  in 
which  Tjc  shows  the  combustion  efficiency. 

CH4  +  202  +  C02  +  2HzO  (16.1) 

CO  +  Io2^C02  (16.2) 

H2+^02^H20  (16.3) 

+  /  Cpjdr)  x%.  =  !>(*+/  CpjdT') 

(16.4) 


3.5.  Modeling  algorithm 

For  the  known  design  parameters,  all  the  above  mentioned  re¬ 
lations  were  applied  for  the  thermodynamic  modeling  of  the 


4.  Exergy  analysis 

The  rate  of  exergy  was  found  in  all  flow  lines  of  the  system  from 
known  thermodynamic  states  and  molar  compositions  by  relative 
chemical  and  physical  equations.  The  rate  of  exergy  destruction  in  a 
component  was  then  computed  from  the  respective  exergy  balance 
equations. 

E^in  =  Ed  +  E^out  +  Wout  (17) 

Or  in  terms  of  fuel  and  product  exergy  streams  and  definition  of 
exergy  efficiency,  we  get 

Efuel  =  ^product  +  ^destructed  (18) 

Vex  =  (19) 

£fuel 

The  relations  for  fuel  and  product  exergy  values  are  listed  in 
Table  1  for  all  system  components. 

5.  Economic  analysis 

Both  thermodynamic  and  economic  aspects  are  important  in 
analysis  and  optimization  of  energy  systems.  Several  methods  were 
proposed  for  the  thermo-economic  assessment  of  energy  systems. 
However,  the  first  and  second  law  of  thermodynamics  in 
conjunction  with  economic  factors  provided  us  a  very  powerful  tool 
for  the  optimization  of  energy  systems.  The  system  total  cost  in  this 
paper  is  one  of  objective  functions  which  should  be  minimized.  The 
total  cost  included  capital  and  maintenance  expenses,  the  system 
operating  cost  (fuel  cost),  as  well  as  the  environment  related  costs 
as  described  in  Section  5.1.  The  investment  cost  functions  are  listed 
in  Table  2  [15,21]  for  major  components  in  terms  of  their  design 
parameters.  In  order  to  convert  these  costs  into  cost  rates  per  unit 
time,  the  following  relations  were  used  [27], 
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Fig.  3.  System  thermal  modeling  flowchart. 


Exergy  relations  for  the  SOFC-MGT  system  components. 


Component 

^product 

£fuel 

AC 

Ei  —Et 

WAC 

FC 

Em  -  Eg 

Wfc 

REC 

Eg  -El 

£6,1 -£7 

FPH 

En  -£io 

£6,2  -  £s 

SOFC 

WsOFC 

(£3-  £n  )  +  (£«-£«) 

CC 

% 

AIP 

Es,2-4l 

£3,1 -£3,2 

GT 

Wgt 

£5  -E6 

HR 

E13  -  En 

£7+£8-£i4 

CK 


CK  x  CRF  x  <p 

N  x  3600 


(20.1) 


1(1+0 
(i  +  on  - 1 


(20.2) 


In  the  above  relations  Ck  shows  the  capital  cost  of  K*  compo¬ 
nent,  CRF  is  the  Capital  Recovery  Factor  which  is  a  function  of 
system  lifetime  (n)  and  the  annual  interest  rate  (i).  <p  shows  the 
maintenance  factor  and  N  is  the  total  working  hours  of  the  system 
per  year.  The  benefit  rate  gained  by  using  the  heat  recovery  heat 
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Table  2 

Cost  functions  for  the  major  components  of  the  SOFC-MGT  system. 


Component 

Cost  function  ($) 

Compressor  (centrifugal) 

Cc  =  91 562(Pc/445)°'67 

Gas  turbine  (radial) 

CGT  =  (— 98.3281n(PGT)  +  1318.5)PGT 

Recuperator 

Crec  =  111.6(mHE)0'95 

SOFC  stack 

Csorc  =  (Atot,stack)(2.96rcell  -  1907) 

Inverter 

C,„v  =  100000(Pceii/500)°'7 

Generator 

CGen  =  60(Pgt  -  Pcf95 

Auxiliary  equipment 

Caux  =  O.ICsofc 

Fleat  recovery  exchanger 

Cex  =  8500  +  405A?  85 

exchanger  is  expressed  in  Eq.  (21)  in  terms  of  difference  between 
the  cost  of  heat  recovery  heat  exchanger  and  the  cost  to  provide  the 
same  required  heat  by  a  separate  boiler  which  includes  investment 
cost,  operational  (fuel)  and  maintenance  costs  of  the  required 
external  boiler. 

^Benefit  =  ^Ext, Boiler  —  OlR, exchanger  (21) 


6.  Environmental  analysis 

Environmental  impact  is  one  of  the  major  concerns  in  analysis 
of  energy  systems  which  is  covered  in  the  present  study  through 
the  optimization  of  the  hybrid  SOFC-MGT  system.  As  a  general 
principle,  an  increase  in  thermal  efficiency  of  a  plant  (which 
provides  a  specific  power  output),  lowers  fuel  consumption  and 
decreases  emission  products  such  as  CO2.  Fuel  cells  are  definitely 
in  a  superior  situation  over  the  conventional  power  plants  in 
terms  of  higher  efficiency  and  lower  emissions.  Experiments 
show  that  the  NOx  emissions  from  a  SOFC-MGT  system  (without 
any  external  fuel  injection  to  the  combustor)  can  be  less  than 
1  ppmv  [28],  In  such  systems,  the  entering  fuel  must  have  the 
lowest  sulfur  content  due  to  sensitivity  of  fuel  cell  parts.  Sweet¬ 
ening  the  inlet  fuel  prevents  the  formation  of  sulfuric  acid  and 
other  similar  pollutants  in  exhaust  gases.  In  this  paper,  a  penalty 
cost  was  associated  to  the  rate  of  C02  emission  produced 
which  was  added  to  the  system  total  cost  during  optimization 
procedure. 


^Ext,  Boiler  —  ^-Cap. Boiler  +  Quel  +  Ol&M  (22) 

The  Payback  period  is  defined  as  the  time  period  required  for  the 
net  income  received  from  selling  the  output  of  a  system  (i.e.  Elec¬ 
tricity  and  recovered  heat)  to  recover  the  initial  investment  costs  of 
the  respective  system  [27],  Regarding  the  time  value  of  money,  the 
worth  of  hybrid  system  capital  investment  as  well  as  the  fuel 
consumption  costs  in  the  pth  year  of  running  can  be  estimated  from 
the  following  equation: 


CInvested,p  =  X/m=0CcaP  m  X  ("*  +  !')P  m  +  ^1-1  Quel  x  N 
x  3600  x  (1  +  i)p~m 


(23) 


7.  Optimization 

7.1.  Objective  junctions 

In  order  to  find  the  optimum  set  of  design  parameters  that 
provide  a  trade-off  between  the  cost  and  performance  of  hybrid 
SOFC-MGT  system  with  minimized  environmental  impacts,  a  multi 
objective  optimization  problem  was  defined.  The  first  objective 
function  to  be  maximized  was  the  total  exergy  efficiency  of  the 
system. 


Obj.Func.I  - 


(25) 


where  i  is  the  annual  interest  rate  and  N  is  the  number  of  running 
hours  per  year. 

In  a  similar  way,  the  net  income  from  selling  the  output  elec¬ 
tricity  and  benefit  gained  by  the  heat  recovery  in  the  pth  year  can  be 
calculated  as  follows: 


The  second  objective  function  to  be  minimized  simultaneously  was 
the  total  cost  of  the  system.  This  objective  function  included  the 
capital,  operation  (fuel)  and  maintenance  costs  as  well  as  the 
environmental  penalty  costs  as  follows: 


Qncome.p  =  £Ll  (Wnet  X  N  X  Celec  +  CBenefit  x  iv)  (1  +  if~m 

(24) 

According  to  the  definition,  the  payback  period  (p)  is  when  these 
two  values  in  Eqs.  (23)  and  (24)  equate.  Then  p  can  be  computed 
using  an  iterative  trial  and  error  method. 


Obj.Func.II  =  Y,  Ccap,o&M  +  Cf  +  Cenv  (26) 

As  the  rate  of  C02  production  is  inversely  proportional  to  the 
system  efficiency  (the  first  objective  function),  the  addition  of  the 
latest  term  in  Eq.  (26)  increases  the  weight  of  thermal  objective 
function  over  the  economic  objective  one  as  the  optimum  system 
tends  to  be  a  more  efficient  system. 


Table  3 

Optimization  constraints  and  range  of  variation  of  design  parameters. 
Constraint  Description 


If 5) < 1203  K 
Tcell  <  1500  K 
ACPR  <  5 
S/C  >1.8 
J?ac,t<0.85 

0.3  <  nwi„  (kg  s-1)  <  0.6 

i  (Am-2)  <4000 

ERec  <  0.9e 

Uf<  0.9 

r14  >  350  K 


Turbine  inlet  temperature  limit 

The  maximum  temperature  value  in  fuel  cell 

The  minimum  value  for  the  steam  to  carbon  ratio 

The  maximum  isentropic  efficiency  for  the  compressor  and  turbir 

Range  of  values  for  Air  inlet  mass  flow 

The  maximum  value  of  cell  current  density 

The  maximum  recuperator  effectiveness  value 

The  minimum  exhaust  gas  temperature 


Range  of  variation  of  design  parameters 
From  To 


1203 

1500 


3  5 


0.7  0.9 
0.7  0.9 
350 
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Table  4 

Genetic  algorithm  optimization  tuning  parameters. 


Parameter  Value 


Population  120 

Maximum  number  of  generations  100 

Probability  of  crossover  (Pc)  0.8 

Probability  of  mutation  (Pm)  0.01 

Selection  process  Tournament 


7.2.  Decision  parameters  and  constraints 

The  design  parameters  (decision  variables)  considered  in  this 
optimization  problem  are  the  air  compressor  pressure  ratio,  the  air 
compressor  isentropic  efficiency,  the  gas  turbine  isentropic  effi¬ 
ciency,  the  inlet  air  mass  flow  rate,  the  cells  fuel  utilization  factor, 
the  cells  current  density,  recuperator  effectiveness  and  the  steam  to 
carbon  ratio  in  reformer.  To  limit  the  results  to  a  range  of  feasible 
points  both  technically  and  economically,  a  set  of  constraints 
(Table  3)  were  applied  to  the  range  of  decision  parameters.  The 
assumed  range  of  change  for  design  parameters  is  also  shown  in 
Table  3. 


Table  5 

Verification  of  modeling  Results. 


Parameter 


Modeling  output  Reference  [4]  Difference  (%) 


Cell  voltage  (V)  0.697  0.7 

SOFC  temperature  (°C)  994.8  990 

Turbine  inlet  temperature  (°C)  896.7  900 

Micro-turbine  power  53.2  53.5 

output  (kW) 

SOFC  power  output  (kW)  218.5  209 

Thermal  efficiency  (%)  61.8  64.9 


0.4 

0.5 

0.3 

0.5 


The  fuel  cost  rate  was  computed  based  on  unit  cost  of  natural 
gas  (0.13  $m  3  or  0.004  $MJ-1)  and  the  environmental  penalty  cost 
was  determined  based  on  the  following  relation: 

Cenv  =  CC02  mCo2  (27) 

where  for  Cc 02  the  penalty  cost  (0.0024  $)/(kgC02)  was  considered 
[29]. 

9.  Results  and  discussion 


7.3.  Optimization  algorithm 

The  above  mentioned  optimization  is  a  multi-objective  problem 
that  due  to  its  high  complexity  and  nonlinearity  and  high  number 
of  decision  parameters,  the  use  of  evolutionary  algorithms  seems  to 
be  the  best  choice  as  a  method  of  optimization.  The  evolutionary 
genetic  algorithms  (GA)  are  sort  of  semi  stochastic  search  methods 
that  imitate  the  laws  of  natural  evolution  for  finding  an  optimal 
solution  to  a  given  problem.  This  method  works  in  the  way  that 
after  formation  of  an  initial  population  of  solutions  as  individuals  as 
different  sets  of  design  parameters,  pairs  of  these  individuals  are 
selected  randomly  and  they  compete  in  terms  of  superior  gens  to 
form  next  generations  that  fulfill  the  survival  of  fittest  in  the  nature. 
With  repetition  of  this  procedure  for  a  number  of  consecutive 
generations,  the  results  approach  to  the  optimum  solution  for  the 
system  [15],  The  multi  objective  GA  in  MATLAB  optimization 
toolbox  has  been  implemented  in  the  present  study.  The  parame¬ 
ters  for  use  with  the  genetic  algorithm  are  listed  in  Table  4  where 
the  stopping  criterion  for  the  algorithm  has  been  set  to  a  specified 
number  of  generations. 


9.1.  Model  verification 

In  order  to  validate  the  modeling  output  results,  the  operating 
parameters  were  compared  with  the  corresponding  data  reported 
in  Ref.  [4]  for  the  same  input  parameters.  Table  5  compares  two 
above  groups  of  data  and  their  corresponding  differences  in  per¬ 
centage.  Results  show  that  the  model  was  capable  of  predicting  the 
thermal  performance  of  the  system  quite  precisely. 

9.2.  Optimization  results 
9.2.1.  Pareto  front 

Two  objective  functions  were  defined  for  the  optimization 
procedure.  Maximizing  the  total  exergy  efficiency  and  minimizing 
the  total  cost  are  obviously  in  conflict  with  each  other  and  there 
is  distribution  of  optimal  points  named  Pareto  front  as  depicted 
in  Fig.  4.  All  single  points  on  Pareto  optimal  front  are  optimum 
non  dominated  solutions.  The  point  at  the  farthest  RHS  shows  the 
highest  efficiency  as  well  as  the  highest  total  cost.  This  point  has 
been  labeled  Point  A  in  Fig.  4  with  62.5%  exergy  efficiency  and 


8.  Case  study 

The  hybrid  SOFC-MGT  system  shown  in  Fig.  1  with  260  kW 
power  and  100  kW  heat  output  (as  80  °C  hot  water)  in  full  load 
operation  was  our  case  study  for  the  system  optimization  in  the 
present  paper.  This  configuration  is  similar  to  the  pressurized 
hybrid  system  pioneered  by  Siemens-Westinghouse  in  2000  [28], 
For  the  economic  modeling  q>  =  1.08  and  N  =  8000  h,  the  interest 
rate  I  =  14%  and  the  unit  cost  of  electricity  ceiec  —  0.06  $1<W_1  h-1 
were  selected  for  use  in  Eqs.  (20)— (24)  [25,26], 

Since  the  issue  of  short  SOFC  lifetime  and  power  degradation 
affects  the  economic  feasibility  of  this  system.  Both  the  SOFC  cost 
and  lifetime  estimations  were  made  with  future  predictions  [15].  In 
this  regard,  40000  h  lifetime  and  0.5%  voltage  degradation  per 
1000  h  of  operation  was  considered  for  the  SOFC  stack  [30],  The 
whole  system  lifetime  (n)  was  considered  equal  to  20  years  while 
the  SOFC  stack  is  replaced  every  five  years  with  a  salvage  value  of  1/ 
3  of  the  purchased  equipment  [31],  while  about  20%  larger  SOFC 
surface  area  was  assumed  to  compensate  for  its  voltage  degrada¬ 
tion  during  the  working  time. 


21 
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Fig.  4.  Pareto  front  for  optimization  results. 
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Fig.  5.  Non-dimensional  Pareto  front  for  optimum  point  selection  using  TOPSIS  deci¬ 
sion  making  method. 


20.71  $h  1  cost  rate.  The  point  at  the  farthest  LHS  shows  the 
lowest  exergy  efficiency  as  well  as  the  lowest  total  cost.  This  point 
has  been  labeled  Point  B  in  Fig.  4  with  56.4%  exergy  efficiency  and 
16.44  $h-1  cost  rate.  Points  A  and  B  designate  a  single  objective 
optimization  with  only  maximum  exergy  efficiency  or  minimum 
Cost  in  mind  respectively.  From  the  behavior  of  the  Pareto  front 
as  can  be  seen  in  Fig.  4  although  the  exergy  efficiency  is 
increasing  from  56.5%  to  about  60%,  the  total  cost  rate  only  in¬ 
creases  slightly.  Obviously  less  fuel  consumption  and  emitted 
pollution  costs  weighs  out  the  associated  increase  in  capital  cost 
of  the  system  in  this  solution  points.  More  rise  in  the  exergy 
efficiency  from  60%  to  higher  values  leads  to  a  drastic  increment 
in  the  total  cost  rate. 

9.2.2.  Optimum  point  selection 

In  order  to  choose  the  best  applicable  solution  on  the  Pareto 
front  both  axes  should  be  non-dimensional  making  it  possible  to 
choose  the  most  optimum  point  based  on  the  distribution  of 
points  and  distances  from  A  and  B  extremes  on  two  sides.  To 
achieve  this  goal,  the  Euclidian  method  has  been  employed  for 
non-dimensionalizing  the  total  cost  rate  [32],  The  Pareto  front 
distribution  of  optimal  points  in  non-dimensional  form  has  been 
plotted  in  Fig.  5.  There  are  various  decision  making  methods  for 
selecting  an  optimum  point  from  the  points  on  the  non- 
dimensional  Pareto  distribution,  in  this  work,  the  TOPSIS 
method  has  been  used  to  select  the  best  optimum  solution  point 
[33],  In  this  method  an  ideal  solution  is  the  imaginary  point 
where  each  objective  has  its  optimum  value  regardless  of  the 
other  objective.  On  the  other  hand,  a  non-ideal  solution  is  the 
point  at  which  both  objectives  have  their  worst  values.  The  two 
ideal  and  non-ideal  solutions  are  shown  on  Fig.  5.  The  best  op¬ 
timum  solution  is  then  selected  as  the  point  form  the  Pareto  front 
which  has  the  shortest  distance  to  the  ideal  and  the  longest 
distance  from  the  non-ideal  solution  points.  The  selection  is  made 
for  the  point  with  minimum  closeness  coefficient  (d_/(d“  +  d+)). 
Where  dr  and  d+  stand  for  the  distances  from  the  ideal  and  non¬ 
ideal  points  on  the  non-dimensional  Pareto  front  respectively 
[33],  Following  this  method,  the  point  Labeled  P  in  Figs.  4  and  5 
has  been  selected  as  the  final  optimum  solution.  Which  corre¬ 
sponds  to  an  exergy  efficiency  of  60.7%  as  well  as  the  total  cost 
rate  of  16.95  $h_1  with  1.95  $h  1  added  for  the  environmental 
emission  externalized  social  penalty  cost.  The  payback  period 
estimation  is  equal  to  about  6.3  years  and  the  output  electrical 


energy  cost  (capital  and  operational  investment)  is  equal  to 
0.057  $kW_1  h 

9.2.3.  Effective  design  parameters 

The  variation  of  design  parameters  for  the  range  of  solutions 
presented  by  the  Pareto  front  is  shown  in  Fig.  6a— h.  From  this  figure 
it  can  be  deduced  that  the  cell  current  density  is  the  most  widely 
distributed  parameter  over  its  allowable  domain.  This  parameter 
also  shows  smooth  changes  versus  to  the  objective  functions.  This 
implies  that  the  current  density  is  an  effective  parameter  in  balance 
between  the  objective  functions.  It  should  be  noted  that  the  current 
density  is  the  parameter  that  determines  the  cell  surface  area  as 
well  as  the  SOFC  stack  cost.  Furthermore  due  to  the  fact  that  the 
exergy  destruction  of  the  stack  (as  the  main  part  of  total  exergy 
destruction)  was  proportional  to  its  current  density,  therefore  the 
exergy  efficiency  of  the  system  was  also  dependent  on  the  cell 
current  density. 

The  system  pressure  ratio  and  air  inlet  mass  flow  are  the  next 
varying  parameters  with  distributed  values  over  its  acceptable 
domain  values.  Other  parameters  are  relatively  concentrated  in 
constant  optimum  values  and  shows  that  they  do  not  have  much 
effect  on  the  conflict  between  the  objective  functions. 

The  values  of  design  parameters  for  the  optimized  case  study  are 
listed  in  Table  6.  In  Table  7  the  dependent  parameters  of  the  opti¬ 
mized  system  are  shown  which  was  selected  from  the  Pareto  front. 
In  addition,  the  temperature,  pressure  and  mass  flow  rates  of  all 
system  points  designated  in  Fig.  1  are  listed  in  Table  8  for  the 
optimized  system.  It  is  noteworthy  that  the  specified  steam  to 
carbon  ratio  in  Table  6  requires  53%  of  the  fuel  mass  flow  to  be 
recycled  at  the  SOFC  anode  as  shown  in  Fig.  2. 

9.2.4.  Exergy  destruction 

The  values  of  exergy  destruction  in  the  system  components  are 
listed  in  Table  9  for  the  optimized  system.  This  table  shows  that  the 
main  sources  of  exergy  destruction  in  the  system  are  the  SOFC, 
combustion  chamber  and  the  air  injection  pipes.  This  is  due  to  the 
fact  that  the  major  portion  of  energy  conversion  takes  place  in  the 
fuel  cells.  Furthermore  due  to  high  irreversible  nature  of  combus¬ 
tion  reaction  and  low  effectiveness  of  air  injection  pipes,  a  large 
temperature  difference  was  present  between  the  entering  air  and 
combustion  products. 

9.3.  Sensitivity  analysis 

In  order  to  achieve  more  insight  into  the  studied  optimization 
problem,  the  sensitivity  of  optimum  values  of  design  parameters 
with  changes  in  fuel  unit  cost  and  capital  investment  were  inves¬ 
tigated.  Fig.  7  depicts  the  variation  of  Pareto  front  for  different 
values  of  fuel  unit  cost,  i.e.  0.004  $Mj_1  as  the  base  value  as  well  as 
0.003  and  0.005  $MJ_1  respectively.  Fig.  7  shows  that  with  higher 
values  of  fuel  cost,  the  Pareto  front  moves  upward  and  to  the  right 
hand  side  (with  higher  exergy  efficiencies).  The  upward  movement 
of  Pareto  front  is  due  to  increasing  the  total  cost  as  one  objective 
function.  Furthermore,  it  seems  that  the  variation  of  objective 
function  is  less  sensitive  to  the  fuel  cost  in  regions  of  higher  exergy 
efficiency. 

The  variation  of  Pareto  front  with  change  in  capital  investment 
is  also  shown  in  Fig.  8.  It  is  clear  that  with  less  expensive  equipment 
the  optimal  system  shifts  to  much  higher  efficiency  to  reduce  fuel 
and  environmental  costs.  With  increase  in  capital  investment,  the 
optimal  values  of  design  parameters  moved  to  the  higher  exergy 
efficiency  values  with  the  lower  values  of  current  density  to 
maintain  the  fuel  cost  as  low  as  possible. 

The  optimum  values  of  design  parameters  with  fuel  cost  and 
investment  cost  change  are  listed  in  Table  10  for  a  selected  point 
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Table  6 

Design  parameters  for  the  optimized  system  selected  from 
Pareto  front. 


Design  parameter 
ACPR 

mail, in  (kg  S_1) 

S/C 

i  (Am  2) 

Uf 


♦  Cf=0.005  $/MJ 

•  Cf=0.004  $/MJ 
▲  Cf=0.003  S/MJ 


2 

16 


14 


X 


Dependent  system  parameters  the  optimized  system  selected  from 
Pareto  front 


12  -I - 1 - T - 1 - 1 - 

0.54  0.56  0.58  0.6  0.62  0.64 


Dependent  parameter 


Total  Exergetic  Efficiency 


Exergetic  efficiency  (%)  60.7 

First  law  efficiency  (%)  63.59 

Overall  CHP  efficiency  (%)  87.8 

Total  Cost  rate  (Sir1)  17.01 

CenvCSh-1)  1.95 

CHP  Benefit  (Sir1)  2.0 

Teen  (K)  1253 

TIT  (K)  1147 


The  values  of  temperature,  pressure  and  mass  flow  rate  for  points  in  Fig.  1  for  the 
optimized  system. 


Point 

m(kgs-') 

P  (bar) 

T(K) 

1 

0.41 

1 

288 

2 

0.41 

4.1 

463 

3 

0.41 

3.97 

857 

4A 

0.0361 

3.77 

1253 

4C 

0.3821 

3.77 

1253 

5 

0.4182 

3.58 

1147 

6 

0.4182 

1.1 

905 

7 

0.3973 

1.04 

534 

8 

0.0199 

1.04 

637 

9 

0.082 

1 

288 

10 

0.082 

10.9 

570 

11 

0.082 

10.57 

794 

12 

0.3 

4 

288 

13 

0.3 

4 

368 

14 

0.4182 

1.01 

381 

from  each  Pareto  front.  Except  for  the  cell  current  density,  the  op- 

timum  values  of  other  parameters  did 

not  show  noticeable  or 

meaningful  variation.  As  explained,  current  density  had  the  most 

prominent  effect  i 

on  the  objective  functions  while  the  pressure  ratio 

and  inlet  mass  flow  had  lower  relative  change. 

9.4.  Change  in  system  size  (power  output) 

In  order  to  study  the  effects  of  change  in  the  system  size  (or 

the  system  power  output)  c 

>n  optimum  design  parameters,  the 

Exergy  destruction  in 

system  components  for  the  optimized  system. 

Component 

Ed  (kW) 

EDIEDitot(%) 

AC 

11.1 

10 

GT 

8.74 

8.2 

SOFC 

36.72 

34 

CC  +  AIP 

44.8 

42.4 

REC  +  FPH 

3.97 

3.6 

Rec.  exchanger 

2.04 

107.37 

100 

Fig.  7.  The  sensitivity  of  Pareto  front  shape  to  fuel  unit  cost. 


optimization  procedure  were  performed  for  200-1200  kW  net 
power  output.  For  each  power  output,  a  proper  optimal  point  was 
selected  and  the  corresponding  results  for  relative  exergy  effi¬ 
ciency  and  system  cost  as  two  objective  functions  are  shown  in 
Fig.  9.  It  was  noticeable  that  with  increasing  the  system  size,  the 
electricity  cost  per  kW  output  reduced.  Another  interesting  point 
was  that  the  optimal  design  parameters  were  very  close  and 
almost  constant  for  different  system  sizes.  The  only  varying  pa¬ 
rameters  were  the  air  compressor  ratio  and  mass  flow  rate. 
Therefore  with  increase  in  the  system  size,  the  optimum  values  of 
fuel  inlet  flow  as  well  as  the  size  of  all  the  system  components 
increased  which  lead  to  a  higher  need  for  oxygen  and  excess 
cooling  air  in  the  fuel  cell.  Furthermore  the  optimum  values  of 
pressure  ratio  increased  for  a  gas  turbine  as  its  flow  rate 
increased.  The  change  of  varying  optimal  parameters  for  different 
sizes  of  hybrid  SOFC-MGT  system  is  shown  in  Figs.  10  and  11. 


10.  Conclusions 

A  hybrid  SOFC-MGT  system  was  modeled  and  optimized  for 
various  power  outputs  (sizes)  using  multi-objective  genetic  algo¬ 
rithm  optimization.  The  objective  functions  were  selected  as  the 
total  exergy  efficiency  and  the  system  total  cost.  To  consider  the 


♦  Capital  Cost=150% 

•  Capital  Cost=100% 
a  Capital  Cost=50% 


12 

10 


•** 


0.56  0.57  0.58  0.59  0.6  0.61  0.62  0.63 

Exergetic  Efficiency 


Fig.  8.  The  sensitivity  of  Pareto  front  shape  to  the  capital 
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Table  10 

Sensitivity  of  the  optimum  values  of  design  parameters  to  fuel  and  capital 


Parameter  C/  =  0.003  C/  =  0.004 


ACPR  4.2  4.1 

t;AC  0.82  0.84 

mairjntkgs^1)  0.42  0.41 

jjct  0.82  0.82 

Uf  0.81  0.82 

f-Yec  0.88  0.89 

i  (Am*2)  2980  2700 

S/C  2.54  2.6 


0.005 


0.83 

0.43 

0.84 

0.79 

0.88 

2560 

2.7 


3.97 

0.82 

0.40 

0.83 

0.83 

0.89 

2445 

2.52 


100% 


0.84 

0.41 

0.82 

0.82 

0.89 

2700 

2.6 


4.21 

0.83 

0.41 

0.81 

0.82 

0.88 

3070 

2.43 


environmental  effects,  a  term  was  added  to  the  total  cost.  It  was 
observed  that  among  design  parameters  of  the  system,  the  cell 
current  density  had  varying  nature  in  the  allowable  range  of  vari¬ 
ation.  This  parameter  played  the  main  role  in  maintaining  a  balance 
between  the  system  cost  and  performance.  However  other  pa¬ 
rameters  such  as  the  pressure  ratio  and  air  inlet  flow  rate  had  next 
level  of  changes  than  that  for  current  density.  The  cell  current 
density  was  also  sensitive  to  the  fuel  unit  cost  and  capital  invest¬ 
ment  changes.  Furthermore  the  change  in  system  power  output 
changed  optimum  values  of  pressure  ratio  and  air  inlet  mass  flow 


Plant  Size  (KW) 


Fig.  10.  The  optimum  values  of  pressure  ratio  for  various  sizes  of  SOFC-MGT  systems. 


rate  significantly  while  the  optimum  values  of  other  parameters 
remained  relatively  constant  with  system  power  output. 
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